Gamma-ray burst: probe of a black hole 
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ABSTRACT 

There is strong evidence for the existence of black holes(BHs) in some 
X-ray binaries and most galatic nuclei, based on different measuremental 
approaches, but black holes aren't finally identified for the lack of very firm 
observational evidence up to now. Because the direct evidence for BHs should 
come from determination of strong gravitational redshift, we hope an object 
can fall into the region near the BH horizon where radiation can be detected. 
Therefore the object must be compact stars as neutron stars(NSs), then 
the intense astrophysical processes will release very high energy radiation 
which is transient, fast- variant. These phenomena may point to gamma-ray 
bursts(GRBs) observed. And recent observations of iron lines suggest that 
afterglows of GRBs show the similar property in active galatic nuclei(AGNs), 
implying GRBs may originate from intense events related to black holes. In 
this letter, a model for GRBs and afterglows is proposed to obtain the range of 
gravitational redshifts(z s ) of GRBs with known cosmological redshifts. Hence, 
we provide a new method that with a search for high-energy emission lines(X- 
or 7-rays) in GRBs, one can determine the gravitational redshift. We expect 
z g > 0.5 or even larger, so that we can rule out the possibility of other compact 
objects such as NSs, suggesting that the central progenitors of GRBs should be 
black holes. 

Subject headings: black hole physics - accretion, accretion disks - gamma-rays : 
burst 



- 3- 



1. Introduction 

The approaches to search for observational evidence for the existence of black holes 
are provided with the various measurement al methods, such as stellar dynamics(Ghez et al. 
1998, 2000), optical emission lines from gas disks (Ferrarese, Ford & Jaffe 1996; Macchetto 
et al. 1997), water maser disks(Miyoshi et al. 1995), X-ray lines(e.g. Fe Ka line, Nandra 
et al. 1997) and the strength of ultrasoft component of X-ray spectra(Zhang, Cui & Chen 
1997). However, all these methods are only related to gaseous processes where(at least two 
gravitational radius) the gravitaional field is weak. To probe the very strong gravitational 
field, the radiation should result from the process near the BH horizon where a large number 
high-energy photons can be released. If the high-energy emission lines are identified, we 
can determine the large gravitaional redshift. Here, we propose that GRBs may be the best 
approach. 

Our method is partly motivated by the recent observational evidence for the existence 
of Fe Ka lines in the X-ray afterglows of GRB 990705 (Amati et al. 2000), 990712 (Frontera 
et al. 2001), 991216(Piro et al.2000) and 000214(Antonelli et al. 2000), and in fact, GRB 
970508(Piro et al. 1999) and 970828 (Yoshida et al. 1999) have also shown the evidence for 
the iron lines. According to the observed line fluence, the mass of iron should account for 
about 10" 4 — 0.1M o , where M Q is the mass of the sun, which contradicts with the standard 
fireball model for GRBs because no enough irons in the interstellar medium can be around 
the progenitors and if the irons are produced during the bursts, it will involve the famous 
problem of bayon contamination. To explain the iron lines, some researches have proposed 
some mechanisms: with the energy injection, only a small mass of Fe is required(Rees 
& Meszaros 2000) and Fe comes from supernova which requires that supernova explosion 
precedes the GRB event by several months to years (Antonelli et al 2001). 

However, the existence of iron lines shows us the some similarity of GRB afterglows 
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to AGNs, implying that the origin of GRBs and their afterglows may be related to black 
holes and accretion process. In this letter, we propose a model of GRBs, in which the 
massive black hole captures a neutron star to produce a GRB and a normal star forming 
an accretion disk to produce the afterglow. Then an alternate mechenism is also proposed 
to explain the iron line origin: Fe line may come from the disk formed by the normal star. 

In the following section, we will give a brief discription of our GRB and afterglow 
model where we will show that some GRBs may be the good candidates of black holes. In 
Section 3, the method to probe the firm evidence for black holes is provided, summary and 
outlook are presented in Section 4. 

2. GRB and Afterglow Model 

Our GRB model is simply described as follows. A massive black hole catches a 
neutron star, which can produce a large number of 7-ray photons through the intense 
astrophysical process in a short timescale. In this case, we could observed only the 
burst. And the minimum variability timescales of burst light curves are related to the NS 
size, At > ^(1 + z)(l + z g ), of the order of milliseconds which is consistent with GRB 
observations, where R is the NS radius, c is light speed and z is the cosmological redshift. 
If a binary (a neutron star and a normal star) rather than an isolated neutron star falls 
into a black hole, when the normal star is disrupted forming a disk and accreted into the 
black hole, X-ray and optical radiation are released. Then both bursts and afterglows 
are observed. Of course, only normal stars captured by black holes in galactic centre will 
produce flares in optical or X-ray band without bursts, which also may be of observability 
(Rees 1988, 1990). And our group now is trying to survey those optical flares. 

The afterglows will be very similar to the continuum radiation of AGNs, for example, 
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the spectrum of the afterglow of GRB 970508 (Galama, T.J. et al. 1998a) may like that of a 
blazar. What's more, the accretion process can account for the afterglows decaying in time 
as a power-law. After the tidal disruption of a normal star occurs, two processes will supply 
mass to the central black hole. The first has been studied by Rees (1988, 1990): the stream 
of stellar mass strung out in far-ranging orbits, showing that the infall rate declines as t~ 5 / 3 . 
The other case involves mass loss from teh inner edge of the accretion disk. Cannizzo et 
al.(1990) have discussed this case of late-time evolution of M in which the accretion disk 
supply rate varies as t~ 12 . 

BH-NS interaction has been studied by some previous authors with the numerical 
simulations (Lattimer & Schramm 1976; Kluzniak & Lee 1998; Janka et al. 1999), which 
was also taken as one of successful models of GRBs. Here, we also have analyzed the 
process with simple calculations. Assuming the different masses of black holes, we derived 
the gravitational redshifts by computing the gravitational radius r g (= ^f^, where G is 
gravitational constant) and the critical radius r where a neutron star is disrupted by the 
central black hole through tidal force. In calculations, we suppose the mass and radius 
of the neutron star are 1.4 M Q and 10 km respectively. For a comparison, the process is 
computed with Landau potential (Landau & Lifshitz 1975): 

0="^ln(l-^) (1) 
and pseudo-Newtonian potential (Paczyhski & Wiita 1980): 

*=--^- (2) 

r-r g 

in separation. In Figure 1, two dashed lines have shown z g as a function of the BH masses 
according to Landau(up) and pseudo-Newtonian(down) potential respectively. We note 
z g ^> 1 when masses are very large, implying that massive black holes may be the better 
test for observations. 
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From the above analysis, a black hole capturing an isolated neutron stars produces 
a large number of high-energy photons through a very intense process, which is a GRB 
observed. With the known total energy of GRBs, we evaluated the gravitational redshift 
around the critical radius 

1 + Zg < -—, (3) 

where E NS is the total rest energy of a neutron star and E y is the isotropic 7-ray energy. 
Since 1 + z g — (1 — ^f)^ 1 ^ 2 , we can find r — r g ~ ( 1+ J 9 - )2 _ 1 - Hence, we estimate the duration 
of a burst, 

T --^ (1 + ^~iow (i + ^-r (4) 

where M is the mass of the black hole, which also tells us the shorter bursts with the 
smaller central masses. The Equation(l) can give a low limit of z g : 

Then the determination of z g will depend on our constraint on the masses of black holes. 
With Eqs. 3 and 5, we can also give a constraint on the masses of black holes, 

— < 10 5 T 7 (-^ - -L). (6) 

When a binary is falling into a black hole, the normal star will be disrupted to form a 
accretion disk in the zone far away from the black hole which is very similar to the accretion 
disks in AGNs, and afterglows in X-ray, optical and radio bands can be produced in the 
disk from where Fe Ka lines are also emitted. We consider a thick disk model, then the 
radial velocity v r ~ ac s , where a is the viscosity parameter (0 < a < 1) and c s is the sound 
speed. The viscosity timescale is given by t v i s ~ ^ ~ where r* is the disk radius. c s can 
be estimated from c 2 s ~ where k is the Boltzmann constant, T is the disk temperature 
and m p is the mass of the proton. Since the radiation of afterglows is emitted from the disk, 
one can find (Frank, King & Raine 1992) T > ( 47r ^a o . ) 1 ^ 4 , where L is the disk luminosity also 
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the afterglow luminosity, and a is the Stefan-Boltzmann constant. Therefore, we reach an 
important result on the masses of black holes: 

where we have taken a ~ 1, and r* ~ 3r g . 

To further test our model, we firstly calculated the isotropic 7-ray energy (£' 7 ) and 
optical (R band) peak luminosity (L) of 20 well-studied GRBs with known cosmological 
redshifts (Bloom, Kulkarni & Djorgovski 2000), and derived some intinsic parameters of 
GRBs which are dispalyed in Table 1. In the calculation, we accept the cosmological model 
with H = 65km s^Mpc -1 ,!^ = 0.3, = 0.7, taking t vis ~ t, the timescale when the 
lightcurves of afterglows evolve from the peak flux to the half. With these parameters and 
above equations, we can estimate the ranges of the black hole masses and z g in our model. 
And R band lumilosity also gives a check on BH masses using the Eddington luminosity 
limit: 

Ml „ k . (8) 

M Q 1.3 x lO^ergs" 1 v ; 

Our final results are presented in Table 1 and plotted in Figure 1, in which two dashed 

lines have been noted above. From Figure 1, we may draw following three conclusions. 

First, we have found gravitational redshifts of some GRBs are very large: seven GRBs 

whose gravitaional redshifts z g > 0.5, which strongely implies that the central engines of 

these GRBs should be the best candidates of black holes. Second, the masses of black holes 

mostly distribute around 1O 6 M which are comparable to the masses of black holes in the 

normal galaxies as in our Galaxy. Third, the range of z g of GRBs calculated is consistent 

with the range given by two lines, supporting our model. 

Because short bursts (with the duration T 7 < 2s) correspnd to the relatively small 
masses of black holes, radiation from the disk around low-mass black holes will be so faint 
that it cannot be observed by our optical telescopes. Therefore, our model for GRBs 
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and afterglows may give a reliable explanation why afterglows cannot be observed from 
most detected GRBs but only in a very small part with the longer durations. Because in 
galaxies the number of isolated neutron stars is larger than the number of neutron stars 
in binary systems according to the observations of pulsars (Taylor, Manchester & Lyne 
1993), the possibility of the former captured by black holes will be much larger, implying 
the afterglows of GRBs are really rare, which is consistent with the observations (Lamb 
2000). The different actions between long-duration GRBs and AGNs may result from the 
circumstance of black holes, since AGNs are in a dirty environment while GRBs in a clean 
one with a neutron star or binary captured occassionally. We think that our unified model 
as a possible mechanism of GRBs and afterglows can successfully explain the different 
occuring rates of bursts and afterglows and the formation of iron lines recently discovered. 

3. Searching for Firm Evidence for Black Holes 

In the above section, we have applied our model to some GRBs with known cosmological 
redshifts. In Table 1, we note that some GRBs have very large gravitational redshifts 
such as GRB 970828 (z g > 3.0), 980425 (z g > 1.4), 000301c (z g > 1.9). Therefore, in our 
model, these GRB central engines should be the good candidates of black holes. However, 
to further probe the central body of a gamma-ray burst, we expect that high spectral 
resolution detectors can find the reliable emission lines in the spectra of GRBs and also 
suggest that pair annihilation lines (511 keV line) may be the best approach. Hitherto, the 
emission lines have been detected and identified in the energy spectrum of GRBs by Venera 
11 and 12 (Mazets et al. 1981) before 1980 and BASTE (Briggs et al. 1997, 1999) in 1990's. 
The centroid energies of emission lines distribute in two different ranges, 330-460 and 40-50 
keV. 

The lines in the broad range of 330-460 keV were interpreted as the strongly redshifted 
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511 keV annihilation lines, but the emission lines at 40-50 keV were thought to be formed 
in another process. For instance, GRB 790526 revealed a emission line at 45 keV which 
was previously thought to be the reliable evidence of emissions at the cycloton frequence 
in the surface of a nearby magnetized neutron star (Mazets et al. 1981). However, since 
the measurement of the redshift of GRB 970508 (Metzger et al. 1997), the evidence 
for a cosmological distance scale for most or all bursts is comfirmed, then GRBs cannot 
be a local phenomina on the surface of a neutron star. Hence, we need a new physical 
mechanism to explain the emission lines observed. In this letter, we interprete the emission 
lines as strongly redshifted 511 keV annihilation lines in the gravitational field of black 
holes. Supposing the cosmological redshift of GRBs, z ~ 1, we calculated the gravitational 
redshifts of the lines whose centroids are around 50 keV, z g ~ 4 in accordance with our 
model expection. Because the gravitational redshift produced by the known compact 
objects including neutron stars or even possible strange stars cannot be over 0.5, we may 
make our decision that the central bodies of these bursts are the best candidates of black 
holes. In fact, the other lines at 330-460 keV observed have the gravitational reshifts in 
the range 0.1-0.5, which is also very difficultly realized on the surface of a neutron star. 
Therefore, GRBs may really originate in the region with the very strong gravitational field, 
where the central progenitors are black holes with the captured neutron stars. 

4. Discussion and Summary 

Because the iron lines discovered in the X-ray afterglows of some bursts imply the 
relation between AGNs and GRBs, a possible model of GRBs and afterglows is proposed 
in the letter. In our model, we can determine the range of gravitational redshifts of 7-ray 
radiation with the intrinsic parameters of GRBs and afterglows, and we note some of them 
are very large (z g > 0.5) implying GRBs as probes of ultra-strong gravity of black holes. 
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However, very firm evidence should come from the identification of GRB line emissions. 
Through the determination of the very strong gravitational redshift of emission lines, one 
can identify the central bodies as black holes. Thus, the spectral line observation of GRBs 
can provide a very powerful evidence for the existence of black holes. Because of the 
level limit on the present missions, we need more advanced detectors with higher spectral 
resolutions in future. Here, we hope that recently launched HETE-II and Swift (in 2003) 
will do much work on the search for GRB emission lines 

The model of GRB afterglows in this letter is also used to explain the X-ray flares 
recently detected in several nearby normal galaxies with the ROSAT data base. Assuming 
the number of massive stars which can produce NSs through supernova explosion is about 
10% of total star number and one supernova event occurs in 100 years, we can simply 
estimate the density ratio of stars and NSs in a normal galaxy as n stax /riNs ~ 1000. The 
observations have given the GRB rate -Rqrb ~ 10~ 6 — 10~ 7 /yr/gal, so we can approximately 
obtain the flare event rate -Rfl are ~ 10~ 3 — 10~ 4 /yr/gal which is consistant with the 
theoretical expected tidal disruption rate (Rees 1988, 1990) and the searching results 
given by ROSAT all-sky survey (Komossa & Dahlem 2001). With the flare event rate and 
assuming a solar mass is swallowed by the black hole per event, we find that the central 
black hole masses in normal galaxies can increase at least up to 10 6 — 10 7 M Q (such as the 
massive black holes in the center of our Galaxy and M31) through the tidal disruption of 
stars within the Hubble time scale (~ 10 10 yr). This conclusion is also important to the 
formation and evolution of normal galaxies and massive black holes. 
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Fig. 1. — The range of GRB gravitational redshifts estimated from our model are plotted 
as a function of the lower limits of the BH masses. Two dashed lines describe the critical 
gravitational redshift of neutron stars disrupted by tidal force according to Landau(up) 
and pseudo-Newtonian(down) potential in separation. 20 GRBs with known cosmological 
redshifts are displayed, and we noticed most masses distribute around 1O 6 M except GRB 
980425 and 990123(as noted in Table 1). If we cancel the optical flash and only take the 
second afterglow of 990123, we note that 990123* which is also plotted in the figure become 
a normal one similar to others. Refer to the text for details. 
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Table 1. The intrinsic parameters of 20 GRBs with known cosmological redshifts and 
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rnnHpl 


name 


T 7 (s) 


logEry (erg) 


logt(s) 


logL(erg s 1 ) 


logM(M Q ) 


z 9 


z 


970228 


47 


51.9 


4.5 


43.2 


5.7-9.2(5.1) 


0.05-358 


0.695 


970508 


18.7 


51.8 


4.5 


45.1 


5.9-8.9(7.0) 


0.16-443 


0.875 


970828 


5 


53.2 


5.0 


45.6 


6.3-6.9(7.4) 


3.0-15.5 


0.935 


971214 


6.8 


53.4 


4.3 


45.8 


5.8-6.9(7.7) 


0.6-10.2 


3.42 


980326 


2.5 


51.5 


4.0 


44.5 


5.4-8.3(6.4) 


0.6-1000 


1 


980329 


13 


54.2 


4.9 


45.6 


6.2-6.3(7.4) 


0.8-1.1 


3.5 


980425 


31 


47.9 


6.0 


45.2 


6.8-13.0(5.1) 


1.4-3X10 6 


0.0085 


980613 


24 


51.8 


4.7 


44.1 


5.9-9.1(6.0) 


0.16-466 


1.096 


980703 


46 


53.1 


4.7 


44.6 


6.0-8.0(6.4) 


0.1-20.5 


0.966 


990123 


15 


54.3 


1.6 


50.2 


4.1-6.0(12) 


0.01-0.3 


1.61 


990123* 






3.5 


46.3 


5.4-6.0(8.0) 


0.15-0.3 




990506 


65 


53.9 


3.7 


45.7 


5.7-7.3(7.6) 


0.05-2.3 


1.310 


990510 


28.6 


53.1 


4.3 


46.6 


5.6-7.5(7.5) 


0.1-19 


1.62 


990705 


36 


52.1 


4.7 


43.0 


5.7-8.9(4.8) 


0.1-214 


0.25 


990712 


21 


51.1 


4.3 


45.7 


5.6-9.6(6.4) 


0.1-2000 


0.434 


991208 


41 


53.2 


4.7 


45.3 


6.0-7.9(7.2) 


0.15-17.7 


0.7055 


991216 


30 


53.9 


4.3 


45.7 


5.7-7.0(7.6) 


0.1-2.5 


1.02 


000131 


9 


54.1 


4.5 


45.3 


5.8-6.2(7.2) 


0.4-1.5 


4.5 


000301c 


3.3 


52.6 


4.5 


45.9 


5.9-7.3(7.8) 


1.9-61 


2.0 


000418 


14 


52.7 


4.9 


44.7 


6.1-7.9(7.6) 


0.6-55 


1.11854 


000926 


8.2 


53.4 


4.3 


46.2 


5.8-7.0(8.1) 


0.4-9 


2.066 



Note. — The timescales T 7 and t have been devided by a factor (1 + z) correction to the GRB intrinsic timcscales. Since a 
particular prompt optical flash (Akerlof et al. 1999) was observed 15 seconds after GRB 990123 stated, we take two cases of the 
afterglow: the flash(m,fl ~ 9) as a peak and the second aftcrglow(mu ~ 17) as the peak where ran is is the visual magnitude 
of optical peak in R band. Because of the possible association between GRB 980425 and SN1998sw (Galama ct al. 1998b), z is 
very small making the 7-ray energy too low, and t is very long probably due to the effect of the light curve of SN1998sw. Mr, 
is shown in the brackets, in many cases, it shows a supcr-Eddington accretion, but some are the sub-Eddington ones. 



